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This Chapter contains a review of the recent results, both experimental 
and theoretical, related to optical control of carriers confined in semi- 
conductor quantum dots. The physics of Rabi oscillations of exciton 
and biexciton occupations, as well as time-domain interference experi- 
ments are discussed. Next, the impact of carrier-phonon interaction in a 
semiconductor structure is described and modern methods of theoretical 
description of the carrier-phonon kinetics and of the resulting dephasing 
are presented. 



1.1. Introduction 

The progress of semiconductor technology that took place in the 80s and 
90s of the last century, in particular the rapid development of epitaxy and 
lithography techniques, has allowed physicists to manufacture and study 
structures in which carriers are confined to a small volume in space (tens 
of nanometers or even less). 1 In this way quantum dots (QDs), that is, 
artificial structures (boxes) containing a few particles with quantized energy 
levels, have been produced. Because of the similarity to natural atoms, such 
structures are also referred to as artificial atoms. However, the properties 
of QDs are more flexible in comparison with atoms: their shapes, size and 
various other features can be engineered at the stage of manufacturing by 
modifying the technological conditions and the number of confined electrons 
can be changed under laboratory conditions (e.g. by applying an external 
voltage) within a wide range of values. 

Parallel to the laboratory investigations, a rapid process of miniaturiza- 
tion of commercial semiconductor structures (e.g., computer chips) took 
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place. At this moment, the state-of-the-art commercially available mi- 
crochips, manufactured using the immersion photolithography technique, 
are built of elements whose size can be reduced down to the 90 nm diffrac- 
tion limit. The introduction of the 45 nm process is announced for 2007 
or early 2008 and the implementation of the 16 nm technology is envis- 
aged in the time frame 2013-140 Thus, the characteristic size of elements 
in the microchips of our standard computer equipment has dropped below 
the size of the first QDs obtained in laboratories 20 years ago and rapidly 
approaches the size of the smallest structures described in today's research 
papers. 

This progress of manufacturing technology is accompanied by a rapid 
development of optical spectroscopy. Currently, it is possible to study the 
optical properties of a single QD and to coherently control the evolution of 
a single carrier or a pair of carriers (electron-hole pair) in such a structure. 
Many experimental schemes of quantum optics have been implemented 
on QDs. Moreover, certain procedures relying on the specific structure 
of the energy levels of these artificial semiconductor structures have been 
demonstrated, which have no counterpart in atomic systems. These experi- 
mental achievements have motivated theoretical proposals for sophisticated 
quantum-optical schemes that may lead, for instance, to optical control of 
a single electron spin in a QD. 

The goal of this chapter is twofold: First, to introduce the reader into the 
fascinating world of modern optical experiments on semiconductor quantum 
dots and to the astonishingly simple, yet nontrivial, theory underlying the 
phenomena observed in the labs on the fundamental, quantum-optical level. 
Second, to give a review of some more sophisticated theoretical methods 
that allow one to include the interaction with the lattice vibration modes 
(phonons) which are specific to semiconductor systems. 

1.2. Essential properties of quantum dots 

Quantum dots are semiconductor structures in which the carrier dynamics 
is restricted in all three dimensions to the length scales of several or a few 
tens of nanometers. 2,3 Various structures that have this property may be 
obtained by a variety of methods. 

One of the most widely used techniques is the Stransky-Krastanov self- 
assembly. When a semiconductor compound is epitaxially (layer by layer) 
grown on a substrate with a different lattice constant (the InAs/GaAs pair 
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is a typical example) each new layer must be squeezed to match the lat- 
tice constant of the substrate. At some point (at about 1.7 monolayers for 
InAs/GaAs) it is energetically favorable for the epitaxial layers to restruc- 
turize into a system of islands, which increases the free surface but relaxes 
strain. 4 The sample is then covered with the substrate material, leading 
to lens-shaped (or, sometimes, pyramidal) nanostructures as in Fig. II. la . 5 
In general, the band edges of the nanostructure material are offset with 
respect to those of the substrate. Here we will only discuss structures as 
those in Fig. Il.lb -d. where the conduction band edge is shifted down and 
the valence band edge is shifted up, so that both electrons and holes are 
bound in the QD structure. 

Another type of structures commonly used in the optical experiments 
are thickness fluctuations of a thin epitaxial layer of a semiconductor (so 
called quantum well) placed between thick structures of different semicon- 
ductor with a wider band gap. Here, again, the band edge offset leads 
to localization of carriers within the quantum well layer. If the epitaxial 
growth of the quantum well layer has been stopped after the formation of 
a new monolayer started, the quantum well has one-monolayer thickness 
fluctuations which weakly localize the carriers. 

The QD nanostructures are typically 2-3 orders of magnitude larger 
than atoms. However, the effective mass of carriers in a semiconductor is 
often considerably lower than the free electron mass (e.g. m* — 0.07mo in 
GaAs), and this degree of confinement is sufficient for quantization of carrier 
energies with electron inter-level spacing reaching 100 meV in self-assembled 
structures. This is definitely enough to resolve the states spectrally and to 
neglect thermal transitions to the excited states even at moderate tempera- 
tures. Therefore, we will restrict the discussion to the ground state of each 
kind of carriers. 

The exact properties of the quantum states in a QD, e.g., the geome- 
try of wave functions or Coulomb interaction between the confined carri- 
ers, may be found, e.g., by tight-binding or pseudopotential calculations. 
For lens-shaped QDs, a simple 2-dimensional harmonic model 2 has been 
shown to be a very good approximation. 6 The states of an interacting few- 
particle system may then be found by numerical configuration-interaction 
techniques. 7,8 However, these details are irrelevant for the discussion on 
the general quantum-optical level. Whenever a specific model is necessary 
(Sees. 11.61 and |1.7|) we will use simple Gaussian wave functions. Also in 
these cases, the results do not depend essentially on this choice. 

Even with the restriction to the lowest orbital states, the structure of 
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Fig. 1.1. (a) Schematic plot of a QD. The xy plane is referred to as the structure plane, 
while z is called growth direction and corresponds to the (approximate) symmetry axis 
of the structure. The arrows show the incidence of light corresponding to the three 
diagrams b-d. (b-d) The diagrams showing the transitions allowed by selection rules in a 
III-V QD, induced by light circularly polarized in the structure plane (b: right-polarized, 
c: left-polarized) and by light linearly polarized along z (d). 



energy levels and allowed optical transitions in a QD becomes quite com- 
plicated if the angular momenta of carriers are taken into account. The va- 
lence band in III-V semiconductors is composed of p-type atomic orbitals 
(orbital angular momentum 1), yielding six quantum states (taking spin 
into account) for each quasi-momentum k. Due to considerable spin-orbit 
coupling the orbital angular momentum and spin are not separate good 
quantum numbers and the valence band states of a bulk crystal must be 
classified by the total angular momentum and its projection on a selected 
axis. Thus, the valence band is composed of three sub-bands corresponding 
to two different representations of the total angular momentum J. Out of 
these, the two states with j — 1/2 form a subband which is split-off by the 
spin-orbit interaction. The other four states with j = 3/2 are degenerate in 
a bulk crystal at k = but this degeneracy is lifted by size quantization and 
strain in a QD structure, with the heavy hole (hh) subband (angular mo- 
mentum projection on the symmetry axis (m = ±3/2) lying above the light 
hole (lh) subband (m = ±1/2) in all known structures. The fundamental 
optical excitation of a QD consists in transferring optically an electron from 
the highest confined state in the valence band (thus leaving a hole) to the 
lowest confined state in the conduction band. The interacting electron-hole 
pair created in this way is referred to as exciton (lh-exciton or hh-exciton, 
depending on the kind of a hole involved). Angular momentum selection 
rules restrict the transitions allowed for a given propagation direction and 
polarization of the light beam, as depicted in Fig. 11.11 For instance, ac- 
cording to the selection rules represented in Fig. II. lb . a right circularly 
polarized (cr + -polarized) laser beam can only create an exciton with total 
momentum +1, referred to as "rr+ exciton", in accordance with the angular 
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Fig. 1.2. (a) The four states of the heavy-hole biexciton system and the optical tran- 
sitions between them. The band diagrams show the particles forming each state. The 
arrows in the valence band represent holes, (b) Schematic representation of the two- 
photon resonance between the ground state and the biexciton state, achieved with a 
properly tuned linearly polarized laser beam. Note that the single exciton transitions 
are forbidden. 



momentum conservation (removing an electron with the angular momen- 
tum m is equivalent to the creation of a hole with the angular momentum 
—m). Similarly, a er_ -polarized (left circularly polarized) beam creates only 
a "<7_ exciton" with the angular momentum — 1. Although there are still 
two transitions allowed for a given circular polarization they can easily be 
distinguished since the lh states are well separated energetically from the 
lowest hh states. 

In appropriately doped structures, QDs in the ground state of the sys- 
tem may be occupied by electrons. An optical excitation in this case corre- 
sponds to a transition between a single electron state and a negative trion 
state, i.e., the state of two electrons and one hole confined in a QD. From 
the Pauli exclusion principle it is clear that this transition is possible only 
if the state which is to be occupied by the photo-created electron is free. 
Hence, in the situation of Fig. lb, a heavy hole trion may be created if 
the dot is initially occupied by a "spin up" (m = +1/2) electron but not 
if the electron in the dot is in the "spin down" (m = —1/2) state. This 
suppression of the optical transition depending on the spin of the electron 
in the QD is referred to as Pauli blocking and has been indeed observed 
experimentally. 9 

After exciting an electron from the heavy hole band with a er + -polarizcd 
laser beam (m = —3/2 — * —1/2), as in Fig. II. lb . it is still possible to 
transfer also the m — +3/2 electron to the m = +1/2 conduction band 
state using cr_ -polarized light (Fig. II. lb ). Thus, if the optical processes 
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are spectrally restricted to heavy holes (as is the case in most experiments) 
there are four optically active states linked by allowed optical transitions 
as shown in Fig. 11.2a . The highest state, with two electrons and two holes 
present in the QD, is treated as composed of two electron-hole pairs and 
is called a biexciton. The Coulomb (dipole-dipole) interaction between the 
two excitons shifts the energy of the biexciton state 10 by the binding energy 
Eb, so that the exciton-biexciton transitions are non-degenerate with the 
ground state-exciton transitions. In this way, all four transitions in the 
diagram can be distinguished either spectrally or by polarizations. 

1.3. Coherent control: experimental state of the art 

The recent progress in ultrafast spectroscopy of semiconductor systems 11 ' 12 
made it possible to control and probe the quantum states of carriers con- 
fined in a QD on femtosecond time scales. In particular, high degree of 
control over carrier occupations in various kinds of QD structures has been 
demonstrated. One kind of an experiment consists in measuring the average 
occupation of the QD after a pulse of fixed length but variable amplitude. 
The QD occupation is defined as the probability of finding an exciton in 
the QD after the laser pulse, calculated as the fraction of cases in which an 
exciton was created over a large number of repetitions of the experiment. 
Within the linear absorption theory, the excitation (the QD occupation) 
grows proportionally to the pulse intensity. Obviously, this growth cannot 
be unlimited. When the occupation of the excited level is sufficiently large 
the spontaneous and induced emission processes suppress further increase 
of the occupation. Here, we are interested in the coherent limit, where the 
evolution is induced by a strong (essentially classical) laser field, inducing 
large occupation changes over time scales much shorter than the sponta- 
neous emission time. Then, the system is driven from the ground to the 
excited state in a coherent way, and than back to the ground state via a 
coherently induced emission process. As a result, in an ideal case, the final 
occupation after a pulse should show sinusoidal oscillations between and 1 
as a function of the square root of pulse intensity (referred to as pulse area 
- see Sec ll.4.2l) . Such oscillations, known as pulse area dependent Rabi 
oscillations, have been indeed observed in a range of experiments on single 
QDs. 13-16 Similar effect can also be observed in ensembles of QDs. 17 

In Fig. 11.3a we show the results of such an experiment, performed with 
a single QD-based photodiode structure. 15,16 ' 18 In this experiment the 
QD is placed in an electric field between a pair of electrodes (Fig. 11.3b ). 
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Fig. 1.3. (a) Pulse area dependent Rabi oscillations, (b) Two-photon Rabi oscilla- 
tions. Solid line: experiment, dashed line: theory (without dephasing), assuming a sech 
pulse shape with a pulse length of 2.3 ps and biexciton binding energy Eg = 2.7 meV. 
Reprinted from Ref. 18, S. Stufler et al, Phys. Rev. B 73, 125304 (2006), ©American 
Physical Society 2006. (c) A schematic explanation of the functionality of the QD pho- 
todiode. 

Each electron-hole pair generated by the optical excitation contributes to 
the photocurrent in the structure. If the laser pulse repetition rate is / 
(typically in kHz-MHz range) and the average QD occupation for a given 
pulse area is n, then the repeated pulsing results in a current / = nef. In 
this way, the Rabi oscillations of the average exciton number are reflected 
in the oscillating form of the photocurrent as a function of the pulse area, 
which provides a means of detection much more efficient than optical ones. 
These oscillations are clearly seen in Fig. 11.3b . although some damping, due 
to dephasing, is also visible. 

Apart from the coherent control of exciton occupations, it has also been 
shown that phase control of carrier states in QDs is possible. A laser 
pulse detuned from the exciton resonance cannot induce real transitions 
and, therefore, does not change occupations. Nonetheless, as shown in 
an optical experiment with interface fluctuation QDs, 19 it can shift the 
energy levels via the AC (optical) Stark effect and affect the phases in 
a quantum superposition of empty dot and single-exciton states. Another 
way to control the phases is to drive the system with a slightly detuned laser 
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pulse, which leads to a combination of occupation and phase evolution 20 
(see Sec. 11.4.21) . In this case, phase effects can be observed in the form of 
Ramsey interference fringes. 21 

It is possible to see in an experiment that the quantum state of an 
electron-hole pair in a QD maintains its phase coherence long after the laser 
pulse has been switched off. To this end, one splits the pulse into two parts. 
If the phase of the quantum state in the QD were random when the second 
pulse arrives the phase of the latter would be irrelevant. On the contrary, 
in the experiment one observes oscillations of the final QD occupation as a 
function of the relative phase shift between the two pulses 14 ' 22,23 (the phase 
is shifted by tuning the delay between the pulses with a sub-femtosecond 
accuracy, that is, by a fraction of the optical oscillation period). Such an 
effect is referred to as time-domain interference, since it may be interpreted 
in terms of interference between the probability amplitudes for exciting the 
QD with the first or with the second pulse. A formal treatment of this class 
of experiments will be given in Sec. 11.4.31 

Apart from the fact that such interference experiments demonstrate co- 
herent phase-sensitive quantum control with an amazingly precise timing, 
they are also of interest from a more general point of view. Obviously, in 
spite of the wave-like behavior manifested by the interference effect, a single 
measurement of the QD occupation always yields either or 1, demonstrat- 
ing the particle-like nature of the exciton. Thus, the time-domain interfer- 
ence experiments on QDs demonstrate quantum complementarity between 
the particle-like nature of an exciton and its ability to show quantum in- 
terference. 2 

Controlling a single quantum degree of freedom (an exciton) is just the 
first step towards large-scale nano-optoelectronic and quantum computing 
applications. The next step towards more complex implementations is to 
include coupling between two or more individual quantum subsystems. The 
simplest experimental realization of such a composite system is a biexciton. 
As explained in Sec. 11.21 due to the Coulomb interaction between the two 
excitons, the excitation energy of an exciton in the presence of the other 
one is different than in its absence. Thus, these two transitions can be 
addressed individually. Together with the polarization dependence of the 
allowed transitions (selection rules, see Sec. I1.2[) , this allows one to excite, 
say, the cr + exciton if and only if the the ct_ exciton is present. Such a con- 
ditional control procedure was indeed performed in an experiment, 25 where 
Rabi oscillations on the exciton-biexciton transition were demonstrated. If 
the two excitons are viewed as quantum bits (with and 1 corresponding 
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to their presence or absence in the QD) then such a conditional excitation 
constitutes an implementation of the controlled-NOT gate which is funda- 
mental for quantum commuting. 26 Such a controlled-NOT gate was indeed 
implemented experimentally in a QD system. 25 In a similar way it is possi- 
ble to coherently manipulate a biexciton system in two coupled QDs (with 
one exciton localized in each dot). 27 

It is interesting to see what happens if the frequency to of a linearly 
polarized laser beam is chosen such that the energy of two photons matches 
the biexciton energy, 2uj = 2E — Eq, while the single-exciton transitions are 
detuned by £b/2 (Fig |1.2b ). A linearly polarized pulse is a superposition 
of two circularly polarized components, so that both exciton transitions 
are allowed by the selection rules but neither of them satisfies the energy 
conservation. Perturbation theory would predict that the occupation of 
the biexciton state should grow proportionally to the square of the pulse 
intensity (that is, to the 4th power of the pulse area), as a result of a two- 
photon absorption process. Experimental results 18 presented in Fig. 11.3b 
indeed show such a behavior for low pulse intensities. However, beyond 
this perturbative regime, a pattern of two-photon pulse area dependent 
Rabi oscillations between the ground and biexciton states develops, which 
become almost periodic for large pulse intensities. The theory of such 
coherent phenomena in the biexciton system will be presented in Sec. 11.51 

1.4. Quantum dot as a two-level system 

Let us now proceed to a theoretical description of the quantum evolution 
of carriers confined in a QD and subject to a laser field. We will start 
with the simplest situation when the QD is driven by a circularly polarized 
laser beam (say, <j + ) tuned to the fundamental heavy hole transition, as 
in Fig. II. la (vertical arrow). Then, only two states are involved in the 
evolution: the ground state (empty dot), denoted by |0), and the cr + exciton 
state, denoted |1) (see also Fig. 11.2a ). Therefore, one effectively deals with 
a very simple two-level system. 

1.4.1. General considerations 

The discussion of such systems is made particularly transparent by intro- 
ducing the concept of the Bloch sphere (actually, a ball). The state (pure 
or mixed) of such a system is represented by a density matrix: a 2 x 2, her- 
mitian, and positive definite operator with unit trace. Any such operator 



February 1, 2008 10:18 World Scientific Review Volume - 9in x 6in chap 



10 P. Machnikowski 




1> 



Fig. 1.4. The Bloch sphere. 

can be written in the form 

p=\{I + n-*l (1.1) 

where n is a real three-dimensional vector with n = \n\ < 1, I is the unit 
operator, and er is the vector of Pauli matrices in the basis |0), |1). Thus, 
the state of a two-level system is represented in a unique way by a unit ball 
in a three-dimensional real space (Fig. II .4[) . Since TrtTi = 0, TrI = 2, and 
(n • er) 2 = n 2 I, one finds Trp 2 = n 2 , so that pure states (represented by 
projectors) correspond to unit vectors n and are mapped to the surface of 
the ball. It is easy to see that the point with spherical coordinates ('d, cp) 
represents the state vector cos(i9/2)|0) + e %v sin(t?/2)|l). 

The electric field of the laser beam is conveniently written as 
£ (t) cos(ojt — <f>), where uj is the frequency of the light, <f> is the phase of 
the pulse and £ (t) is the envelope of the pulse, varying slowly compared 
to the oscillations of the optical field. The state of the system \^) evolves 
according to the Schrodinger equation 

ihf t m=H\*), (i.2) 

with the Hamiltonian 

H = E\l){\\ + £{t - i a ) cosMt - Q - 4>] (m|0X1| + M*|1)(0|) , (1.3) 

where E is the energy of the interband transition in the QD, /i is the off- 
diagonal (interband) matrix element of the electric dipole moment (between 
the two relevant states), and we allow the pulse to arrive at an arbitrary 
time i a - By a proper choice of the phase of the state 1 1) it is possible to have 
/i real, which will be assumed in the following. We will use the shorthand 
notation f(t) — fi£(t). 
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In the absence of the driving, the phase of the state |1) rotates with 
the frequency E/H which is of the order of fs _1 . One gets rid of this trivial 
fast dynamics by describing the system in the frame rotating with the same 
frequency. This is similar to a transition to the interaction picture except 
that, for practical reasons, it is convenient to perform the transformation 
with the laser frequency lj instead of the system frequency E/h (these two 
frequencies are close to each other). Thus, we perform the transformation 

|#)=C/ r |*), U t = ex P M|l)(l|). 

Using Eq. ()1.2|) one easily finds the evolution equation for the redefined 
states in the form ih\^>) = H\^f), with the Hamiltonian 

H = -tkj\l)(l\ + U v HUf (1.4) 



-A|1X1| + ~f(i-i. 



e — iuit — i<f> _|_ g—iui(2t — t )+i4> 



e iujt Q +icj) , e iu(2t—t )—i(t> j H 



where A = fko — E is the detuning of the laser beam from the transition 
energy. 

The next essential step is to note that the natural frequency scale of 
the system evolution is set by the detuning A and by the pulse amplitude 
f(t) and is many orders of magnitude smaller than the optical frequency lj. 
Therefore, the quickly oscillating terms (~ e 2luJt ) are strongly off-resonant 
and will have very little impact on the system evolution. Therefore, they 
can be neglected, which leads to the rotating wave approximation (RWA). 28 
As a result, one obtains the following RWA Hamiltonian 



RWA 



= -A|l)(l| + ~/(t- t a ) (e-^ to +«|0)(l| + e^+^lXOl) (1.5) 



= -A|lXl| + ~/(f-ta)«-ff, 

where u — [cos(tjt a + <fi), sin(wt a + 0), 0]. 



1.4.2. Pulse area dependent Rabi oscillations 

Let us start with applying the formalism to the Rabi oscillations described 
in Sec. 11.31 First, assume that the pulse is resonant with the optical tran- 
sition in the QD, that is, A = 0. In this case, the evolution operator 
generated by the RWA Hamiltonian (|1.5[) can be found analytically. In- 
deed, in the resonant case the Hamiltonians at different times commute 
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Fig. 1.5. (a) Pulse-area dependent Rabi oscillations at the resonance. The points 
represent final system states for different pulse areas a £ (0, 2n) with a step of Aa = 0.1. 
(b) Pulse-area dependent Rabi oscillations off resonance. The points represent final 
system states for different pulse areas a £ (0, Git) with a step of Aa = 0.2. (c) The 
occupation of the state |1) as a function of the pulse area for resonant (solid) and non- 
resonant (dashed) driving. 



with one another, so that the evolution operator may be written as 



U(t) = cxp 



2 v ' 



• • *(*) . 
i sin — - — u ■ 



(1.6) 



where = f^drffr), t is the initial time of the evolution, and the 
last identity is easily proven by expanding the exponent in a series and 
collecting the odd- and even-order terms using the identities (u ■ cr) 2m = I 
and (u ■ er) 2m+1 — u ■ cr for any natural number m. We will always assume 
that this time is before any pulses were switched on, so that one can set 
t{) — > — oo. The value of a = $(oo) is called the pulse area and determines 
the unitary transformation of the system state performed by the complete 
pulse. One speaks of 7r-pulses, ir/2 pulses etc., referring to the value of the 
pulse area. 

Starting from the ground system state |0) one obtains, after switching 
the pulse off, 



|*> = tf(oo)|0>=coB-|0> 



■ id) • 

ie v sin 



a 



|1>- 



(1.7) 



The occupation of the state |1) is, therefore, 

Kl|*)| 2 = |<l|l-)| 2 =sin 2 



and indeed oscillates between and 1 as a function of the pulse area a. 
The final system states for a set of values of a in this resonant case are 
shown in Fig. 11.5b . Let us note here that the final state is determined by a 
simple function of just one quantity, the pulse area, and is independent of 
any details of the pulse shape. This fact is known as the area theorem. 
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Off resonance (for A ^ 0), the evolution can be found in a closed analyt- 
ical form only for rectangular pulse envelopes f(t). Instead of performing 
this simple exercise, let us look at the evolution for a family of Gaussian 
pulses 

ha i ( t V 
/(i) = _^ e Hnr) (1.8) 



27TT p 

obtained by (also very simple) numerical integration of the Schrodinger 
equation (| 1 . 2[> . shown in Fig. 11.5b (see also Ref. 20). In the detuned case, 
the system does not reach the |1) state. Instead, the state is rotated around 
a tilted axis and the final state does not show periodicity as a function of 
the pulse area. The difference between the resonant and non-resonant case 
is also visible in the dependence of the final occupation of the state |1) 
(Fig. Ob). 

1.4.3. Time- domain interference 

Another experiment that can easily be explained based on the two level 
model is that of time domain interference. As discussed in Sec. 11.31 such 
experiments are performed with two laser pulses selectively tuned to the 
exact resonance with one of the two fundamental optical transitions. The 
pulses are generated by splitting a single laser pulse and delaying one part 
with respect to the other. In this way the two pulses are phase- locked, i.e., 
their relative phase is definite and determined by the delay time r. The 
first pulse arrives at t a = and prepares the initial superposition state. 
In a usual two-slit (space-domain) experiment, this would correspond to 
splitting the particle path. The second pulse arrives at t a = t. This pulse 
plays the role of "beam merger" providing, at the same time, a phase shift 
between the "paths" . 

The Hamiltonian describing the system driven by the two pulses is ob- 
tained by an obvious generalization of Eq. (|1.5p . 



H = |/a(*)(|0><lj + |1)(0|) + \h{t r)(e— |0)<l|^ r |l)<0|), (1.9) 

where fi{t) are the envelopes of the two pulses and we have set tfi — 0. The 
pulses do not overlap in time so that the evolution can be split into two 
independent stages. 

In an experiment, the system is initially in the state |0). The first pulse 
is a 7r/2 pulse that performs the transformation U\ = (I — ia x )/^/2. This 
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pulse leaves the system in the equal superposition state 



M - M d.10) 

The second pulse is again a 7r/2 pulse, 

U 2 = ^=(I-in-a), (1.11) 
v2 

where n = [cos wr, sin wx, 0] , as follows from the general discussion in 
Sec. 11.4.11 After this pulse, the average number of excitons in the dot 
is 

7V(0) = |(l|f/ 2 |^)| 2 = i(l-cos W r), 

and changes periodically between iV m j n = and iV max = 1 as a function of 
the delay time r, thus producing an interference pattern. 

The quality of the interference pattern is quantified in terms of visibility 
of interference fringes 

N max N m j n 



V 



The amplitude (visibility) of the fringes in an ideal experiment is V = 1 for 
7r/2 pulses, i.e., for an equal superposition state in between the pulses. Oth- 
erwise, some a priori information on the superposition state can be inferred 
and the visibility is reduced. Reduction of visibility occurs also if some in- 
formation on the system state between the pulses (analogous to which way 
information) has been extracted either intentionally, in a controlled exper- 
iment (see Sec J1.5.2| . or as a result of uncontrolled dephasing. 24 

1.5. Beyond the two levels: optically driven evolution in a 
biexciton system 

The quantum control of a two-level system formed by the ground state and 
a single exciton exploits only small part of the possibilities offered by a QD. 
Much more interesting physics, as well as potential applications, becomes 
available if all the states are involved in the dynamics. Below we study two 
quantum optical schemes involving the confined biexciton system. First, 
the theory of the experimentally observed two-photon Rabi oscillations is 
given. 18 Then, a theoretical proposal for an experiment highlighting the role 
of quantum complementarity in time-domain interference is described. 24 
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1.5.1. Two-photon Rabi oscillations 

If we apply circularly polarized excitation to a QD in the ground state we 
can only create single excitons, as follows from the selection rules discussed 
in Sec. 11.21 and from the Pauli exclusion principle. By using linear polariza- 
tion, which is a superposition of two equal cr + and er_ polarized components, 
one enables also biexciton generation. In this case, the ground state is cou- 
pled to the biexciton state via both exciton states. Here we will consider 
a situation in which the laser frequency is tuned to the two-photon res- 
onance with the biexciton, while the singlc-cxciton states are detuned by 
A = -Eb/2, as shown in Fig. 11.2b . 

In the rotating wave approximation, the Hamiltonian of the biexcitonic 
system described above can be written in the form 

H = ^ (l+X+l + I-X-D + 0= [(|g) + |XX))«+| + (-|) + H.c] , (1.12) 

where |+),|— ) and |XX) are the a + , er_ and biexciton states, respectively. 
Thus, the laser pulse couples only one (bright) combination of the states 
ground and biexciton states, |B) = (|g) + |XX))/\/2 to one (-^-polarized) 
single exciton state |X) = (|+) + |— ))/v / 2- 

The system evolution is easily described in the new basis \Y) = (|+) — 
|-»/V2,|D) = (|g) - |XX»/V2, |X),|B), where the first two states are 
the (decoupled) ^-polarized exciton state and the dark combination of the 
states ground and biexciton states. In this basis, the Hamiltonian (|1.12p 
can be written as 



H = 



/E B /2 \ 


E B /2 f(t)/V2 

\ 0f(t)/V2 / 



(1.13) 



which corresponds to the detuned rotation in the invariant subspace 
spanned by the x-polarized exciton state \X) and the state |B), while the 
first two states are decoupled and undergo only a trivial evolution. The 
initial (ground) state is now written as |g) = ^ D ^ B ^ . 

Let us consider the instantaneous (for fixed t) eigenstates of the Hamil- 
tonian (|1.13[) as a function of f(i). For a large binding energy E-q the two 
branches belonging to the non-trivial two-dimensional invariant subspace 
are always widely separated. Thus, under the action of a sufficiently slowly 
varying laser pulse the system undergoes an adiabatic evolution with the 
pulse envelope fit) playing the role of a slowly varying parameter of the 
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Hamiltonian. Therefore, one can assume 30 that at each time t the state 
corresponds to the adiabatic (instantaneous) eigenstate of the Hamiltonian 
(Q51) . i.e., 

= |D) + [c-W |X)+c + (t)|B)] e -" A W 



where 



and 



c ± {t) = -= 1± 



V2 
Eb 



1/2 



A(t) = i ^ rfr (i? B - V(^) 2 + 8/ 2 (r)) . (1.14) 

After the laser pulse has been switched off, the state becomes 
m _ P>±£^H _ ,_„_,„ („^l |g) + „,„ AM |XX) 

(1.15) 

and coincides with the biexciton state for A(oo) = tt. Thus, the system 
effectively performs Rabi oscillations between the ground state and the 
biexciton state, with the occupation of the biexciton state given by 



2 A(oo) 



ATxx = |(XX^)r = sin / ^-A (1.16) 

For 9 «C T AE B /h one has 

» r s 1 f°° , ,2/ n 4?iArcoshv / 2 /1 2 
Moo) « — — / dt/ 2 (t) = 



ftA.Bs 7_ oc ir 2 AE B T 

and the biexciton occupation A^xx grows as # 4 , i.e., proportional to the 
square of the pulse intensity, as expected for a two-photon process (see 
Fig.OJ. 

It should be noted that the pulse area appearing as the parameter of the 
standard Rabi oscillations is now replaced by the adiabatic dynamical phase 
A, which is a nontrivial functional of the pulse shape [Eq. (|1.14|) ]. In this 
way, the simple universal dynamics off a resonantly driven two level system, 
described by the area theorem, is replaced by a more complicated pattern 
of oscillations, depending on the detuning parameter t p Eb but also on the 
exact pulse shape (e.g. Gaussian vs. sech 18 ). The absence of any single 
cxciton occupation in the final state is obviously due to the fact that the 
adiabatic limit corresponds to weak and long (i.e., spectrally narrow) pulses 
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and the excitation of single-exciton states becomes forbidden by energy 
conservation. Nonetheless, during the evolution the single-exciton states 
are also occupied. 

1.5.2. Quantum complementarity in time-domain interfer- 
ence experiments 

This Section discusses the essential modification to time domain interfer- 
ence experiments (Sec. 11.4.31) that allows one to attain partial information 
on the state of the system and to observe the related visibility reduction 
of the interference pattern. First, however, a measure of the partial infor- 
mation is introduced and the complementarity principle is formulated in a 
quantitative form. 31,32 

The notion of "partial information" is understood as follows. The sys- 
tem (S) of interest is coupled to another quantum probe (QP) system and 
conditional dynamics of the latter is induced, leading to correlations be- 
tween the states of the systems S and QP. Next, a measurement on QP is 
performed and its result is used to infer the state of S, i.e., to predict the 
outcome of a subsequent measurement on S. The probability of a correct 
prediction ranges from 1/2 (guessing at random in absence of any correla- 
tions) to 1 (knowing for sure, when the systems are maximally entangled). 

Quantitatively, an intrinsic measure of information on the system S 
extracted by QP is provided by the distinguishability of states, ' 



where p is the probability a correct prediction for the state of S maximized 
over all possible measurements on QP. In this way, guessing at random and 
knowing for sure correspond to T> — and T> = 1, respectively. According to 
a general theory, 31 ' 32 the complementarity relation between the knowledge 
of the system state and the visibility of the fringes may be written, using 
the distinguishability I? as a measure of information, in the quantitative 
form, 



The equality holds for systems in pure states. 

In a QD, this formal scheme translates naturally into the conditional 
dynamics of a biexcitonic system (Fig |1.2a ). as described in Sec. 11.31 Thus, 
the exciton addressed in the interference experiment described in Sec. 11.4.31 




(1.17) 



V 2 + V 2 < 1. 



(1.18) 
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Fig. 1.6. (a) The diagram of energy levels and transitions in the system, assuming 
that the excitons are confined in neighboring dots with different transition energies. The 
pulses inducing transitions on the system S have broad spectrum so that both transitions 
are possible. The pulse acting on the QP system is spectrally selective and tuned to the 
biexcitonic transition, so that the single-exciton transition is energetically forbidden in 
this subsystem, (b) The sequence of pulses used in the experiment. 

is, say, the er+ exciton (system S). The other degree of freedom of the 
biexciton system (the er_ exciton) will be used as the quantum probe. 

We will use a tensor product notation with |0) and |1) denoting the 
absence and presence of the respective exciton, as previously, with the in- 
terfering system (S) always to the left. In the rotating basis with respect 
to both subsystems, the RWA Hamiltonian for the biexciton system is 

H = H 1 ®1+ A|l)(l| <g> |1)(1| + i/Qp(t - T QP )I <8> (|0)(1| + H.c.) (1.19) 

where A is the bi-exciton energy shift and £ , Qp(t) is the envelope of the 
pulse coupled to the QP exciton. Here the first term denotes the Hamilto- 
nian (|1.9[) and corresponds to the pulse sequence of the interference exper- 
iment described in Sec. 11.4.31 the second one accounts for the bi-excitonic 
energy shift and the third term describes the action of the pulse coupled 
to the second (QP) exciton and spectrally tuned to the exciton-biexciton 
transition. This pulse arrives at t = tqp, between the other two pulses 
(that is, < iq>p < t), and will induce the conditional dynamics. Its phase 
is irrelevant and will be assumed 0. The structure of system excitations 
and the sequence of pulses are shown in Fig. 11.61 

Assume that the exciton (system S) is in the equal superposition state 
(jl.lOp . The probability of correctly guessing the result of a measurement in 
the |0), |1) basis without any additional information is obviously 1/2. Now, 
we can correlate this excitonic system with the other one (QP; initially in 
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the state |0)). To this end, one applies a selective (spectrally narrow) pulse 
with the area a (the arrow labelled 'QP' in Fig. II. 6h ). This rotates the 
state of the second subsystem according to Eq. (|1.6[) (with h = [1, 0, 0]) if 
and only if the first system is in the state |1). The corresponding unitary 
transformation of the compound system is 

Uqp = |0)(0| ® I + |1)(1| ® (cos -I i sin -a x ) , 
and takes the state into 

W) = ^10) ® 10) - <8> (cos||0) - ism||l>) . 

For a = 7r, this pulse performs a CNOT-like transformation 26 on the biex- 
citonic system. As a result, the total system is in the maximally entangled 
state (|0)|0) — i|l)|l))/2 and a measurement on the QP system uniquely 
determines the state of the system S. Hence, due to quantum correlations 
between the systems, complete information on the state of S has been ex- 
tracted to QP. On the other hand, if the biexciton is excited with a pulse 
with area a < tt the correlation between the subsystems is weaker and a 
measurement on QP cannot fully determine the state of S, although the 
attained information may increase the probability for correctly predicting 
the result of a subsequent measurement on S. According to the discussion 
above, this means that partial information on the state of S is available. 

In order to find the distinguishability measure in the biexciton scheme 
discussed here, we write the density matrix of the total system correspond- 
ing to the state \ip'), 



\^2\n){m\® Pnm , (1.20) 



B 2 

nm 

with poo = |0)(0|, pn = |(I + cosa u z — sina a y ), p 01 = p{ = 
icos f |0)(0| — sin^|0)(l|. Note that poo and pn (but not p i) ar e density 
matrices. 

According to the general theory, 31 ' 32 the best chance for correctly guess- 
ing the state of S results from the measurement of the observable poo — pn 
and the probability of the correct prediction is then p = | + jTr|poo ~ Pn\ 
where \p\ is the modulus of the operator p. Using the explicit forms of the 
density matrices PoOjPii and the definition (| 1 . 1 T[) one finds for the distin- 
guishability in our case 



v = ^Tr|p o-pii| 



sm- 



(1.21) 
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Thus, the amount of information on the system S accessible via a measure- 
ment on QP increases from (no QP pulse at all) to 1 (for a it pulse). 

Next, we study the effect of extracting the which path information on 
the interference fringes. In the state (|1.20jl . the reduced density matrix of 
the subsystem S is ps = TrQpg = (1/2) (I — cos(a/2)a y ). Upon applying 
the second pulse of the interference experiment scheme, namely the uncon- 
ditional 7r/2 pulse [Eq. the average number of excitons in the dot 
is 



Now, the average occupation oscillates between the limiting values (1 ± 
| cosa|)/2 and the visibility of the fringes is 



Comparing Eq. (|1.22[) with Eq. (| 1 . 2 1 1) it is clear that the more certain 
one is whether the exciton is there or not (T> increases), the less clear 
the interference fringes become (V decreases). Quantitatively, the relation 
V 2 + V 2 = 1 holds which is consistent with the complementarity relation 
(|1.18|) . It can be shown that, in the presence of coupling to the environment, 
where both subsystems are in mixed states due to dephasing, this relation 
will turn into inequality. 24 Let us notice that the impact of the which 
path information on interference fringes is the same no matter whether the 
QP subsystem is measured before or after generating and detecting the 
interference fringes and even whether it is measured at all. 

The time-domain manifestation of quantum complementarity discussed 
here not only broadens the class of experiments in which fundamental as- 
pects of the quantum world may be tested but also has the advantage of 
being independent of the position-momentum (Heisenberg's) uncertainty 
that has been historically tied to the space-domain discussions of com- 
plementarity. 33 In fact, it is independent of any uncertainty principles 
whatsoever. Indeed, the only two quantities which are measured in the 
time-domain interference experiment are the occupations of two different 
excitons. These quantities refer to different subsystems and, therefore, are 
obviously commuting and simultaneously measurable. Although the quan- 
tum probe exciton is created in a way that correlates it with the presence 
of the other exciton (S), this cannot be interpreted as a (projective) mea- 
surement on S, since the QP exciton is definitely a quantum (microscopic) 
system and not a classical, macroscopic measurement device. Thus, the 
experimental procedure described in this section demonstrates quantum 




V = |cos(a/2)|. 



(1.22) 
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complementarity in its pure form, involving only the notion of information 
on the system state and independent of any uncertainty relations between 
non-commuting observables. 

Further discussion of the concept of complementarity in the context of 
optical experiments on QDs, including the analysis of its feasibility in terms 
of the currently available experimental techniques and of the parameters of 
the existing structures, as well as the analysis of the role of dephasing is 
given in Ref. 24. 

1.6. Carrier-phonon interaction in quantum dots 

This section presents the derivation of the coupling constants between bulk 
phonon modes and confined carriers in a semiconductor. 34,35 We will re- 
strict ourselves to the deformation potential coupling, which is relevant for 
the discussion in the following sections. The treatment of other couplings 
(piezoelectric and polar) is reviewed in Ref. 36. 

Any crystal deformation leads to shifts of the conduction (c) and va- 
lence (v) bands which are, to the leading order, proportional to the relative 
volume change. The corresponding contribution to the energy of electrons 
(e) and holes (h) in the long-wavelength limit is 

where <7 e /h are the deformation potential constants for electrons and holes 
and V is the unit cell volume. Using the strain tensor cr, 

1 / dv.i duj 

one may write H^^' — ^fa c / h Tra — T^c/h^ ■ u(r), where u(r) is the local 
displacement field. The displacement is quantized in terms of phonons, 

^)=^E/5^^(^+^)^. (^) 

where Wk is the frequency for the wave vector k, = —e~k is the corre- 
sponding real unit polarization vector, and p is the crystal density. Only the 
longitudinal branch contributes to V • u in p.23[) and the final interaction 
Hamiltonian in the coordinate representation for carriers is 
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In the second quantization representation with respect to the carrier 
states this reads 



nn' 

i 



n\V{r)\n')a\a n > = 



nn' 



d 3 ril)*(r)V(r)iJj n >(r)ala n > 



^ knn' 



(1.25) 



where f^ nn ,{k) = <r c/h sj ^h^nn* (k), with the formfactor 
,(fc) = / d 3 rV>;(r)e ife - r Vv(r) - 



(1.26) 



While the common coefficient of the coupling Hamiltonian contains the 
fundamental and material-dependent constants and reflects general electri- 
cal and mechanical properties of the semiconductor system, the formfactor 
(|1.26j) contains information about the geometry of the confinement and the 
resulting properties of wave functions. In this sense, it is the "engineerable" 
part of the carrier-phonon coupling. 

From the orthogonality of single-particle states one immediately has 
•7"nn'(0) = 8 nn i. If the wave functions are localized at a length I, then the 
extent of the formfactor is ~ 1/7. Thus, for carrier states localized in a QD 
over many lattice sites and smooth within this range, the functions T nn >{k) 
will be localized in fc-space very close to the center of the Brillouin zone. 

As an example, let us consider a Gaussian wave function, 



1 



7T 3 / 4 U_L 



■ exp 



(1.27) 



where r± is the position component in the xy plane and l± , l z are the local- 
ization widths in-plane and in the growth (z) direction. The corresponding 
formfactor is then easily found to be 



T(k) = exp 



kill 



(1.28) 



The formulas derived above describe the interaction in the single- 
particle basis. However, most of the following deals with excitonic states, 
i.e. states of confined electron-hole pairs interacting by Coulomb poten- 
tials. Both carriers forming the exciton couple to phonons according to 
Eq. (|1.25p . For the further discussion in this chapter it is sufficient to dis- 
cuss the lowest exciton state. It is reasonable to assume its wave function 
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approximately in a product form, i.e., 




11, u 

the QD. Then, one gets from Eq. (|1.25|) the following coupling constant 



between the confined exciton and phonons for the deformation potential 
interaction 



Note that, due to different deformation potential constants Cc/h, this cou- 
pling does not vanish even if the electron and hole wave functions are the 
same, leading to identical single-particle formfactors. 

Thus, the exciton (restricted to its ground state) driven by a laser field 
and interacting with LA phonons via the deformation potential coupling is 
described by the Hamiltonian 



where the first term describes the carrier subsystem and is given by 
Eq. (HI). 

The Hamiltonian ()1.30|) is the basis for the microscopic modeling of 
dephasing effects in QDs. It turns out that this relatively simple model is 
very successful in reproducing the experimental data. 38,39 Hence, it may 
serve as a reliable starting point for describing the evolution of the combined 
system of confined carriers and lattice modes. 

1.7. Theoretical methods for carrier phonon kinetics 

In this section we will study a few theoretical methods that have proven 
to be useful for the description of the carrier-phonon quantum kinetics in 
QDs. This discussion will be limited to the simplest case of a two-level 
system described by the Hamiltonian (|1 .30|) . 

1.7.1. Exact solution for ultrafast excitations 

Let us start with the case of an ultrafast excitation. A very short laser 
pulse prepares the system in a certain superposition (dependent on the 
pulse phase and intensity) of the ground state |0) (no exciton) and the 
single-exciton state |1). By very short we mean a pulse much shorter than 




(1.29) 



H = i?R WA + ^kbibk + 1 1 >< 1 1 9k(b k + bl k ), 



(1.30) 



k k 
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the time scales of phonon dynamics, so that the preparation of the initial 
state may be considered instantaneous. This corresponds to the actual 
experimental situation with pulse durations of the order of 100 fs. 40 On the 
other hand, the pulse is long enough to assure a relatively narrow spectrum 
and to prevent the population of higher confined levels and excitation of 
optical phonons. 41 In this ultrafast limit the only role of the laser pulse is 
to prepare the initial state, while the subsequent evolution is generated by 
a time-independent Hamiltonian, 

H = -A|1}(1| + J2 *» h blb k + |1X1|$>*(&* + ftlfc), (1-31) 

k k 

obtained by setting H RWA = -A|l)(l| in Eq. (fl~30|) . 

In a superposition state created by the laser pulse [Eq. Q 1.7)1 ] the inter- 
band component of the electric dipole moment has a non- vanishing average 
value oscillating at an optical frequency (hence referred to as optical polar- 
ization) 42 which leads to the emission of coherent electromagnetic radiation 
with an amplitude proportional to the oscillating dipole moment. In an un- 
perturbed system (e.g., in an atom), the radiation would be emitted over 
time of the order of the lifetime of the superposition state, i.e., until the 
system relaxes to the ground state due to radiative energy loss. 

In a semiconductor structure an additional effect, related to carrier- 
phonon coupling, appears on a time scale much shorter than the lifetime of 
the state. The last two terms in Eq. (|1.3ip describe a set of harmonic os- 
cillators which, in the presence of the exciton, are displaced by an external 
force proportional to the coupling constant gk- This means that, due to the 
interactions between confined carriers and lattice ions, the ground state of 
the lattice in the presence of a charge distribution is different than in its ab- 
sence. As a result, after the creation of a confined exciton the lattice relaxes 
to a new equilibrium, which is accompanied by the emission of phonon wave 
packets 37,43 that form a trace in the macroscopic crystal distinguishing the 
exciton state from an empty dot. This information broadcast via emit- 
ted phonons leads to a decay of the coherence of the superposition state 44 
although the average occupations of the system states remain unaffected 
(hence the process is referred to as pure dephasing). Since coherent dipole 
radiation requires well-defined phase relations between the components of 
a quantum superposition, the amplitude of this radiation, measured in the 
experiment, gives access to the coherence properties of the quantum state 
of confined carriers itself. The dephasing of the quantum superposition is 
therefore directly translated into the decay of coherent optical radiation 
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from the system. 

As discussed in Sec. II. 6[ the carrier-phonon interaction term in 
Eq. (|1.31| is linear in phonon operators and describes a shift of the lat- 
tice equilibrium induced by the presence of a charge distribution in the 
dot. The stationary state of the system corresponds to the exciton and the 
surrounding coherent cloud of phonons representing the lattice distortion 
to the new equilibrium. The transformation that creates the coherent cloud 
is the shift wbkW^ = bk — fk/(huJk), generated by the Weyl operator 35,44 



A straightforward calculation shows that the Hamiltonian (|1 .31|) is diago- 
nalized by the unitary transformation W = |0)(0| (8)1+ |1)(1| ®ro, where I is 
the identity operator and the tensor product refers to the carrier subsystem 
(first component) and its phonon environment (second component). As a 
result one gets 



We assume that at the beginning (t = 0) the state of the whole system 
is p = ( | \E r )(\E f | ) ® p p h, where p p h is the density matrix of the phonon sub- 
system (environment) at thermal equilibrium and |^) is given by Eq. (jl.7p . 
For simplicity, we will assume an equal superposition state, setting a = n/2 
and tp = — 7r/2. 

The evolution operator U (£) = e -' lHt / h m ay be written as 

U(t) = W^WU(t)W^W = W^U(t)W = wW{t)wu\t)u{t) 
= W^W(t)U(t), 

where U(t) = e - iEt l h and Wit) = U(t)WW(t). Since U(t) is diagonal the 
explicit form of W(t) may easily be found and one gets 



where w(t) — e lH ph t / h we lHpht ^ h . 

Using the evolution operator in the form (|1.33|) the system state at a 
time t may be written as 




(1.32) 



H = WHW^ = -A|l)(l| + H ph , 
where A = A + J^k \9k\ 2 /{hw k ) and H ph = J2k ^kb\b k . 



U(t) = [|0)(0| ® I + |1)(1| ® w^w{t)} U(t) 



(1.33) 




(1.34) 
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where we used the tensor product notation in which an operator A is ex- 
panded as A = J2 m n \ m )( n \ ® A mn with a set of operators A mn acting on 
the second subsystem, and written in the matrix form with respect to the 
first subsystem. The density matrix for the carrier subsystem is obtained 
by tracing out the phonon degrees of freedom, i.e., ps = TrEp. Hence, 



Ps(t) 



1 

2 y - 



1 e- iAt / h (w^(t)w) 
iAt / h {w^w(t)) 1 



(1.35) 



The average may be calculated using two rules for Weyl operators. 34,44 The 
multiplication of two operators of the form 



exp 



yields 



W1W2 — W3 exp 



5 Eft" 



• (2) 



i = 1,2,3 



(1) (2)« 

Tk % 



(1.36) 



(1.37) 



where W3 is given by Eq. (|1.36[) with 7^ 3 -* = 7^ + 'y^ ■ The rule for 
averaging of an operator given by Eq. (|1.36[) in the thermal equlibrium 
state is 

( u; .) =e -iEj7< <) | 2 (2n fc +l) ; (138) 

where n\~ are bosonic equilibrium occupation numbers. The final result is 



exp ■ 



E 



fjk 



[i sinwfct + (1 — coscjfet)(2nfc + 1)] 



The emitted coherent dipole radiation is proportional to the non- 
diagonal element of the density matrix p s (t) and its amplitude is 



P{t) = P \(wHt)w)\ 



(1.39) 



In Fig. 11.71 we show the normalized polarization amplitude P(t)/Pg (orig- 
inally derived in Ref. 45). The interaction with the macroscopic crystal 
environment leads to a reduction of coherent radiation due to pure dephas- 
ing of the exciton state, reflected by the reduced value of the non-diagonal 
element of the density matrix ps. At t = one has (w'(t)w) = 1, while 
at large values of t, coswfct oscillates very quickly as a function of k and 
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Fig. 1.7. Decay of the coherent radiation from a confined exciton at various temper- 
atures, as shown. In our calculations we use typical parameters for a self-assembled 
InAs/GaAs structure: single particle wave functions ip(r) modelled by Gaussians with 
= 4 nm and l z = 1 nm [Eq. Ql.27E ]. the deformation potential difference rj e — cr^ = 9.5 
eV, crystal density g = 5300 kg/m 3 , and the speed of sound c = 5150 m/s. 

averages to (see also Ref. 45). Thus, for long times, the polarization 
amplitude tends to a temperature-dependent finite value 

2 



This partial decay of coherence is a characteristic feature of short-time 
dephasing for carrier-phonon couplings encountered in real systems. 

1.7.2. Perturbation theory for driven systems 

In this section we derive the equations for the reduced density matrix of the 
carrier subsystem to the leading order in the phonon coupling, assuming 
that the unperturbed evolution is known. This approach yields very simple 
and intuitive formulas that may easily be applied to a range of problems. 
Here, we will again restrict the formalism to a single exciton level. A more 
general discussion is given in Ref. 36. 

As already mentioned, the evolution of the carrier subsystem is gener- 
ated by the (time dependent) Hamiltonian -Hrwa [Eq. (11. 5(1 ], describing the 
properties of the system itself as well as its coupling to driving fields. The 
evolution of the phonon subsystem (reservoir) is described by the Hamil- 
tonian ffph [the second term in Eq. (|1.30|) ]. The evolution operator for 
the driven carrier subsystem and free phonon modes, without the carrier- 
phonon interaction, is 
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where U^f\{t) is the operator for the unperturbed evolution of the carrier 
subsystem and s is the initial time. 

The carrier-phonon coupling may be written in the form V — S ® R, 
where S acts in the Hilbert space of the carrier subsystem while the time- 
independent R affect only the environment. For instance, in the special 
case of Eq. (TOO)) . S = |1)(1| and 

R =^T,^{b k + bl k ). (1.40) 



N k 

We will assume that at the initial time s the system is in the product 
state g(s) = IV'oXV'ol <8> /0 p h, where \ipo) is a certain state of the carrier 
subsystem and p p h is the thermal equilibrium distribution of phonon modes. 
Physically, such an assumption is usually reasonable due to the existence 
of two distinct time scales: the long one for the carrier decoherence (e.g. 
1 ns ground state exciton lifetime 40 ' 46 ) and the short one for the reservoir 
relaxation (1 ps dressing time 37,40 ' 45 ). 

The starting point is the evolution equation for the density matrix of 
the total system in the interaction picture with respect to the externally 
driven evolution Uq, in the second order Born approximation with respect 
to the carrier-phonon interaction 47 

Q(t) = + ^f* dr[V(r), g(s)\ - ± £ dr jT dr'[V(r), [V(r'), g(s)]}, 

(1.41) 

where g(t) = U^(t)g(t)U (t), V(t) = U%(t)VU (t) (it should be kept in 
mind that, in general, V may depend on time itself). 

The reduced density matrix of the carrier subsystem at time t is p(t) = 
^RWa(*)p(£)^rvva(*)' P(t) = TrRg(i), where the trace is taken over the 
phonon degrees of freedom. The first (zeroth order) term in (|1.4ip gives 

P {0 \t) = cwWIMV-oKwaW = \Mt))(Mt)l (i.42) 

The second term vanishes, since it contains the thermal average of an odd 
number of phonon operators. The third (second order) term describes the 
leading phonon correction to the dynamics of the carrier subsystem, 

P (2) (t) = -i fdr [ T dr'TT R [V(r),[V(r'),e( S )}}. (1.43) 



& .. 

First of the four terms resulting from expanding the commutators in 
G33J is (I) = -Qt|MV>o|, where 

Qt = ¥j s dT J s dT ' s ( T ) s ( T ')( R ( T - T ') R )- I 1 - 44 ) 
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The operators S and R are transformed into the interaction picture in the 
usual way 5(f) = U$(t)SU (t),R(t) = U&(t)RU (t) and (6) = Tr R [Op ph ] 
denotes the thermal average (obviously [Uo(t), p p h] = 0). 
The second term is 

(II) =^JdrJ dr'|^o)(^o|^(r')5(T)(i?(r' - t)R) = -\^o\Ql 

where we used the relation (R(t' — t)R) = (R(t — t')R)* . 

In a similar manner, the two other terms may be combined into 

(m) + (iv) = $t OM^oO- 

where 

h\p] = Jpf dr J dT'S(T')pS(T)(R(r-r')R). (1.45) 
In terms of the new Hermitian operators 

M = Qt + Qt, ht = ^(Q t -Qt), (1-46) 

the density matrix at the final time t (jl.42ll.43[) may be written as 

P(t) = U RWA (t) (\ip )(M ~ * [ht, \ipo)(M} (1-47) 



-\ {A t , iMV'ol} + *t[\MW]) ul WA . 

The first term is a Hamiltonian correction which does not lead to irre- 
versible effects and, in principle, may be compensated by an appropriate 
modification of the control Hamiltonian Hrwa- The other two terms de- 
scribe processes of entangling the system with the reservoir, leading to the 
loss of coherence of the carrier state. 

Let us introduce the spectral density of the reservoir, 

R(") = / dt{R(f)R)e iwt - (1-48) 
For the operators given in Eq. (| 1 .40|) one has explicitly 

R(w) = in B (w) + Ifffel 2 P(w - wik) + S(uj + w fc )] , (1.49) 

fc 

where ub(w) = — u>) — 1 is the Bosc distribution function. The spec- 
tral density R(u>) depends on the material parameters and system geometry 
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and characterizes the properties of the lattice subsystem. For the deforma- 
tion potential coupling to LA phonons, it has the long- wavelength behavior 
R(w) ~ Lu 3 [n B {uj) + 1]. 

With the help of (|1.48j) one may write 



$ t [p] = J dojR{uj)Y(uj)pY^(uj) 
where the frequency-dependent operators have been introduced, 

Y(uj) = J drS{T)e luJT . 
Using (|1.48p again one has 

Qt = J duj j dr j dT'0{ T - T l )S{T)S{T l )R{uj)e- tuj{T - T 
Next, representing the Heaviside function as 

6{t) = -e u 

we write 



(1.50) 



(1.51) 



^—itj't 



dioR(uS) 
= - I du)R(w) 



2ni uj' -uj + i0+ 
dm' Y^{u')Y{oj') 

-inS(ui' - w) + V — 



2iri to' -uj + i0+ 
duj 1 
2ni 



Y\u')Y{J) 



where V denotes the principal value. 

Hence, the two Hermitian operators defined in (|1.46p take the form 



Ai 



and 



duR{uj)V 



<koR{uo)Y\Lo)Y{u) 

du'Yi(w')Y(u)') 



(1.52) 



(1.53) 



2tt uj 1 — uj 

In order to provide an example of an application of the perturbative 
method presented in this section let us consider a 7r/2 rotation of an ex- 
citonic two-level system in the presence of carrier-phonon coupling. As 
we have already pointed out, any fast change of the state of the carrier 
subsystem leads to spontaneous processes of lattice relaxation that affect 
the coherence of the carrier state. It is reasonable to expect that coherent 
control is recovered if the evolution of the carrier subsystem is slow (adi- 
abatic) compared to the typical timescales of the lattice dynamics. Thus, 
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the requirement to avoid traces of carrier dynamics in the outside world 
favors slow operation on the carrier subsystem, contrary to other decoher- 
ence processes (of Markovian character) , like radiative decay of the exciton 
or thermally activated processes of phonon-assisted transitions to higher 
states. The latter have the character of an exponential decay and, for short 
times, contribute an error 6 — T g /rd, where r g is the gating time and is 
the decay time constant. Here we consider the interplay between these two 
contributions to the error for the solid-state qubit implementation using ex- 
citonic (charge) states in quantum dots (QDs), 48 with computational states 
defined by the absence (|0)) or presence (|1)) of one exciton in the ground 
state of the dot, operated by resonant coupling to laser light. We show that 
it leads to a trade-off situation with a specific gating time corresponding to 
the minimum decoherence for a given operation. 

To quantify the quality of the rotation, we use the fidelity 2 ® 



which is a measure of the overlap between the ideal (pure) final state with- 
out perturbation, L/r,wa(oo)V'0j and the actual final state of the system 
given by the density matrix p(oo). If the procedure is performed ideally, 
i.e. without discrepancies from the desired qubit operation, then F = 1. 
The fidelity loss 6 = 1 — F 2 is referred to as the error of the quantum gate. 
Using the definition (|1.54|) and the Master equation (|1.47p . the error may 
be written in a general case as 



It should be noted that the unitary correction generated by h t does not 
contribute to the error at this order. Using the definitions (|1.45)1.52|) this 
may be further transformed into 



where ip± is a state orthogonal to in the two-dimensional space of 
interest. 

Since the coherence of superpositions induced by short pulses is unsta- 
ble due to phonon- induced pure dephasing (Sec. I1.7.T]) . it seems reasonable 
to perform operations on dressed states, i.e. on the correctly defined quasi- 
particles of the interacting carrier-phonon system. 50 This may be formally 
achieved by employing the solid-state-theory concept of adiabatic switch- 
ing on/off the interaction 51 (as done in Ref. 52) to transform the states of 



F = (■ip \Ul WA (oo)p(oo)U R wA(oo)\ip Q ) 1 



(1.54) 





(1.55) 
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the nonintcracting system into the states of the interacting one. Thus, we 
assume adiabatic switching on/off of the interaction with phonons by ap- 
pending the appropriate exponent to the original interaction Hamiltonian 
[Eq. CE3DM, 



ii>(ii5>fc ( b " +b 



fct 

-fc 



(1.56) 



k 

where e = + . The operator S now becomes 

S(t) = U RWA (t)e-^\l)(l\Ul WA (t), 

where the free evolution operator is generated by an optical pulse at the 
resonance (Sec. 11.4. 2[) . The general formula (|1.55p may now be used with 
the bare initial state \ipo)- The adiabatic procedure assures that it is trans- 
formed to the dressed state before comparing it to the density matrix p, so 
that the fidelity is defined with respect to stable, dressed states. 
The operator Y(oj) can be written in the form 

Y(w) = ^FHdlXOl - |0)(1| + |0)(0| - |1)<1|) 

+ ^F*(- W )(|0)<1|-|1)(0| + |0)(0|-|1)<1|), 



where 



/CO J 
-co « T 



Since in quantum information processing applications the initial state 
of the quantum bit is in general not known, it is reasonable to consider the 
error averaged over all input states. Let us introduce the function 



S(u)=oj 2 \(iP ± \Y(uj)\tPo)\ 



2 

av ' 



where the average is taken over the Bloch sphere. According to Eq. (j 1 . 55[) , 
the average error may now be written as 

5 = J ^R{cj)S(lj). (1.57) 

The averaging is most conveniently performed by noting that 

Y{UJ) = i F(w)l+>H + i F * ( ^ )hX+l ' 
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where |±) = (|0)±|1))/V2. Choosing \%) = cos § |+)+e^ sin § |-), = 
sin 1 1 +) — e lv cos 1 1 — ) , one gets 



5 H = \ 



F(-co)e ttp cos 2 - - F*(uj)e- iv sin 2 - 



2 



which, upon averaging over the angles 9, if on the Bloch sphere, leads to 
S(cj) = ±(\F(uj)\ 2 + \F(-u>)\ 2 ). 

Let us now consider a Gaussian pulse for performing the quantum gate, 
f(t) = [a/(v / 27TTp)]e _ 2 (*/ r p) ) were r p is the gate duration, while a is the 
angle determining the rotation. The function |F(w)| 2 that carries all the 
needed information about spectral properties of the system's dynamics may 
be approximately written as 



|^±H| 2 ~a 2 e Tp V V2 "W . (1.58) 

As may be seen from (|1.5T|) and (|1.58[> . for a spectral density R{tu) ~ ui n 
the error scales with the gate duration as t~" +1 and t~" +2 at low and 
high temperatures, respectively. Therefore, for n > 2 (typical for various 
types of phonon reservoirs) the error grows for faster gates. Assuming the 
spectral density of the form R{to) = Rdp^ 3 for low frequencies (as for the 
deformation potential coupling at low temperatures), we obtain from (|1.57|) 
and IL581) 

5 = -^a 2 R DP T- 2 , atT = 

This leading order formula holds for S <C 1. Also, if we introduce the upper 
cut-off, the error will be finite even for an infinitely fast gate (see Fig. 11.81) ; 
this is the ultrafast limit discussed in Sec. 11.7.11 

Thus, the phonon-induced error indeed decreases for slow driving. This 
could result in obtaining arbitrarily low error by choosing a suitably low 
gate speed. However, if the system is also subject to other types of noise 
this becomes impossible. 49 Indeed, assuming an additional contribution 
growing with rate 7m, the total error per gate is 



7nM . 1 2 n 1 , i n> \ 

-^2-+7mt p , 7 nM = —a i? D p, 7m = — , (1.59) 



where r r is the characteristic time of Markovian decoherence (recombination 
time in the excitonic case). As a result, the overall error is unavoidable and 
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Fig. 1.8. (a) Combined Markovian and non-Markovian error for a a = tt/2 rotation on 
a qubit implemented as a confined exciton in a InAs/GaAs quantum dot, for T = (solid 
lines) and T = 10 K (dashed lines), for two dot sizes (dot height is 20% of its diameter). 
The Markovian decoherence times are inferred from experimental data. 40 (b) Spectral 
density of the phonon reservoir R(ui) at these two temperatures and the gate profile S(u>) 
for a = tt/2. 



optimization is needed. The formulas (|1.59|) lead to the optimal values of 
the form (for T = 0) 

3 ( 2giR BP \ 1/3 (2 \ 1/3 

£>min = - ^ J , for T p = \^-a RupTr J . (1.60) 

For the specific material parameters of GaAs, the optimal gate time and 
minimal decoherence resulting from Eqs. (|1.60p are 

t p = a 2 ' 3 1.47 ps, r5 min = a 2 / 3 0.0035. 

The exact solution within the proposed model, taking into account the cut- 
off and anisotropy (flat shape) of the dot and allowing finite temperatures, 
is shown in Fig. 11.81 The size-dependent cut-off is reflected by a shift of 
the optimal parameters for the two dot sizes: larger dots allow faster gates 
and lead to lower error. 

It should be noted that these optimal times are longer than the limits 
imposed by level separation. 48,53 ' 54 Thus, the non-Markovian reservoir ef- 
fects (dressing) seem to be the essential limitation to the gate speed. On 
the other hand, in the above discussion we used simple Gaussian pulses and 
a straightforward way of encoding the qubit. The error resulting from the 
phonon dynamics may be reduced by optimizing the shape of the control 
pulse 55 ' 56 or by encoding the qubit into a state of an array of QDs. 57 
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1.7.3. Correlation expansion 

The correlation expansion technique is a standard method used for the de- 
scription of quantum kinetics of interacting carriers and phonons in semi- 
conductor systems of any dimensionality. 58 It has been successfully ap- 
plied to carrier-phonon kinetics in quantum dots driven by an optical field, 
beyond the instantaneous excitation limit described in Sec. 11.7.11 and be- 
yond the weak perturbation case which allows a perturbative treatment 
(Sec. 11.7. 2p . 59-61 Compared to higher-dimensional systems, in quantum 
dots coherent and non-equilibrium phonons play a larger role because of 
the localized polaron effect. Therefore, a reliable description of the car- 
rier phonon-kinetics in these systems requires a high enough degree of the 
correlation-expansion technique. 61 

Various implementations of this technique differ in notation and in the 
choice of dynamical variables. Here, let us start from the three dynam- 
ical variables x,y,z describing the carrier state, x — {<j x (t)), . . ., where 
<7i(t) — e lHt / h aie lHt / h are Pauli operators, written in the |0), |1) basis, in 
the Heisenberg picture. These three variables are the coordinates of the 
evolving Bloch vector, uniquely determining the reduced density matrix of 
the carrier subsystem according to Eq. (jl.ip . 

From the Heisenberg equations of motion one finds the dynamical equa- 
tions for these three variables, 

x = i{[H,a x ]) = -Ay-4y^ReS fe -4y^Rey fe , (1.61) 

k k 

and analogous for y and z (from now on, the time dependence will not 
be written explicitly). Obviously, this set of equations is not closed, but 
involves new the phonon variables B k = g k (b k ), as well as phonon-assisted 
variables of the form y k = g k ((<J y b k )) = (vybk) — (vy)(b k ). The double an- 
gular brackets, ((. . .)), denote the correlated part of a product of operators, 
obtained by subtracting all possible factorizations of the product. 

Next, one writes down the equations of motion for the new variables 
that appeared in the previous step, for instance, 

ilk = i([H,y k ]) = Ax k - 2Vz k - iuj k y k + \g k \ 2 (iyz + x) (1-62) 
+2 ^2(x qk + x qk ) + 4x k Re B q + 2x^{B qk + B qk ), 

<? 1 Q 

where the new two-phonon and two-phonon-assisted variables arc defined 
as B qk = g q g k ((b q b k }}, B qk = g* q g k {{b q b k )) , x qk = g q g k {{<r x b q b k )), x qk = 
9 q 9k{{&xb q b k )) , etc. In the next step, one writes the equation of motion 
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Fig. 1.9. Rabi oscillations of the exciton occupation /=(l+*)/2asa function of the 
nominal pulse area a, calculated using the correlation expansion at T = 0. Reprinted 
from Ref. 61, A. Krugel et al, Phys. Rev. B. 73, 035302 (2006), ©American Physical 
Society 2006. 



for these new variables, introducing three-phonon variables. It is clear that 
the resulting hierarchy of equations in infinite and has to be truncated at a 
certain level. Here we do this by setting all the correlated parts of the three- 
phonon and three-phonon assisted variables equal to zero. This amounts to 
neglecting the correlations involving three or more phonons or, physically, 
to neglecting three-phonon processes (that is, emission or absorption of 
three or more phonons within the memory time of the phonon reservoir, 
which is of order of 1 ps). The motivation for this procedure is that higher 
order correlations should play a decreasing role in the dynamics. From 
the equations of motion it is also clear that such higher order correlations 
develop in higher orders with respect to the coupling constants gk- 

As an example of an application of this technique, Fig. 11.91 shows the 
results for the Rabi oscillations of a coherently driven exciton (sec Sec. II. 4|) 
interacting with phonons. 60 An interesting feature of this result is that the 
quality of oscillations decreases for moderate pulse durations but then in- 
creases again for longer pulses. This can be understood (using the perturba- 
tive approach of Sec. 11.7.2]) in terms of a resonance between the oscillations 
of the charge distribution and the phonon modes. 62 



1.8. Conclusion 



This chapter presented an overview of some recent results related to optical 
control and decoherence of carriers in semiconductor quantum dots. The 
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problems discussed here are of interest not only from the scientific point of 
view but are also important for possible applications of nanostructure-based 
devices in the field of nano-electronics, optoelectronics and spintronics. 63,64 

Obviously, a chapter of limited length cannot give an exhaustive review 
of this broad and rapidly developing field. Consequently, the goal of the 
present review was rather to give an introduction to the field rather than 
a complete, encyclopedic account of all the achievements. In particular, 
the two level (or few-level) model on which the presented discussion was 
based, is obviously merely an approximation to the complex semiconductor 
system. In fact, any interaction of the confined carriers with the external 
driving fields not only leads to the desired quantum transitions but also can 
induce some unwanted ones. If both the desired and unwanted transitions 
have a discrete nature (e.g., the exciton vs. biexciton transition in a QD) 
the latter may be suppressed by a suitable choice of control pulses. 53 ' 54 A 
more demanding problem of transitions involving the macroscopic phonon 
continuum is treated in Ref. 56 

Another class of experiments and theories that are not covered by this 
review concerns optical methods for coherent spin control in QDs. The 
capability of encoding and manipulating information at the single-spin level 
is of great importance for spintronic and quantum information processing 
applications. Recent experimental progress includes the generation 65 and 
optical control 66 ' 67 of the spin coherence together with a possible read- 
out of the state of a single confined spin in a QD system. 68 It was also 
demonstrated that spin states in QDs may be prepared with high fidelity 
(exceeding 99.8%) by optical coupling of electronic spin states. This was 
done by resonant excitation of the trion transition in the presence of small 
heavy-light hole mixing. 69 Corresponding to these achievements, a range 
of schemes for optical spin control has been proposed theoretically. 70-75 
These schemes are based on quantum-optical procedures, exploiting the 
structure of selection rules in a QD discussed in this chapter. In particular, 
the spin rotation is performed by coupling af a single electron state to a 
trion state. Since the spin control is achieved by means of spin-dependent 
charge evolution these procedures result in phonon-induccd dephasing that 
may be described using the methods described here. An introductory 
review of optical spin control in QDs is given elsewhere. 77 
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